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Abstract 


A propeller lifting surface design and analysis program is automatically coupled with a Euler 
/ Integrated Boundary Layer Theory (IBLT) axisymmetric flow solver. The lifting surface 
method solves the localized propeller problem, while the Euler/IBLT solver handles the global 
flowfield, capturing the effective inflow problem. For viscous flows, the boundary layer is 
constructed based upon the parameters specified by the IBLT solution, and is merged with the 
inviscid Euler flowfield. The robust coupled method is capable of representing open propellers, 
ducted proulsors, and internal flow cases, including multi-blade row applications. For large 
axisymmetic bodies, the user may specify a nominal inflow, and the coupled method is used 
for the localized propulsor problem only, further increasing the computational efficiency. The 
specified nominal flow field may be calculated by other numerical flow solvers, obtained from 
experimental results, or calculated from a Euler/IBLT solution of the entire body. The coupled 
code is an extremly rapid flowfield gridding, calculation, and convergence method, which allows 
an order of magnitude reduction in convergence time when compared to the current. efforts 
using Reynolds Averaged Navier Stokes (RANS) as the flow solver. Experimental validation 
is demonstrated for open, ducted, and internal flow propulsors. 
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Chapter 1 


Introduction 


1.1 Overview 


The goal of this research was the creation of a fully automatic coupled method between a 
lifting surface propeller blade design and analysis program and an Euler throughflow solver. 
A coupled technique is necessary as the propeller influences the flowfield, and the flowfield 
influence the propeller. Thus, both the propeller and flow field calculation must be iteratively 
converged to a final solution. Currently, coupled design and analysis of marine propulsors is 
generally accomplished by using a Reynolds Average Navier - Stokes (RANS) flow solver, the 
use of which is a computational intensive task requiring a high degree of operator knowledge 
and experience. A Euler solver with an Integral Boundary Layer Theory (IBLT) boundary 
layer reconstruction routine, enables solution convergence twenty to thirty times as fast as the 


current RANS methods. 


1.2 Lifting-Surface Propeller Blade Design and Analysis 


The use of lifting surface theory for the design and analysis of propellers has been extensively 
documented in the literature, and will be only summarized here for the sake of completeness(3},[11],[12. 
Propeller Blade Design (PBD-14) developed by Kerwin! et al was used for this research. A 


lattice of discrete vortex segments is placed along the blade mean camber surface, the hub, the 


‘Justin E. Kerwin, Profesor of Naval Architecture, Ocean Engineering Department, MIT. 
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duct, if present, and the trailing wake system. Control points are placed a particular point 
within each grid lattice. Each lattice segment has an appropriate vorticity strength, such that 
the kinematic boundary conditions are satisfied along the gridded sivace, namely the total 
velocity must be tangent to the surface at each control point. 


In propeller design and analysis, the velocity relation is 


=> = => 
4 total = V inducea + V ef fective (1.1) 
where 
Vee = the total velocity in the presence of an operating propeller 
Vacs = the velocity induced by the propeller and trailing wake vortex distributions 
rae ffective = the velocity due to the interaction between the propeller vortex distributions 
and the presense of vorticity in the inflow 
and 


Ve coal = the velocity present in the absence of the propeller 


=> => 
It is important to realize that neither V inducea nor V ef fective are physical velocities that 


may be measured. Only Vis actually exits in the flow field, and its division into two 
components is done solely for computational reasons. It is important to point out that the 
above listed velocities are completely independent from the nominal velocity Ves which 
is the velocity present in the absence of the propeller. Based on the above definitions, it is 
useful to note that if no vorticity is present in the nominal flowfield, then no interaction between 
the propeller, as modeled by vortex distributions, and the inflow are necessary. Thus, in this 
special case, the nominal inflow is equal to the effective inflow. This is never the case in reality 
since viscous effects are always present, however it often serves as a useful check of propeller 
analysis code operation when run in an inviscid mode. 


The lifting surface method solves the propeller problem by first developing an influence 


13 


matrix [INF] which gives the velocity induced by each vortex lattice segment, with an assumed 
unit strength, upon each control point. When this influence matrix is multiplied by the actual 


vortex segment strength, the induced velocity at every control point is known. 


[INF] [Tr] = Vlas (1.2) 


The kinematic boundary condition is specified by dictating that the total velocity normal 


to the grid surface must be zero. Thus, for a zero thickness case, 


[INF] [r| ae V es fective| -A=0 (1.3) 


A lifting surface code may be used for either a design of a new propeller or an analysis of 
an existing propeller. In both cases, Equation 1.3 applies, only the knowns versus unknowns 
vary. In the case of propeller design, a desired radial and chordwise loading distribution is 
prescribed, and the blade shape is manipulated until Equation 1.3 is satisfied. In propeller 
analysis, the blade shape is prescribed, and the resulting circulation is solved. Blade thickness 
effects are added by placing discrete source lines coincident with the spanwise blade vortex 
lines. The resulting propeller forces due to the lifting surface and thickness are calculated 
from the Kutta-Joukowski and Lagally theorems, respectively. A leading edge suction force 
and Lighthill pressure distribution correction are applied to those forces, and the propeller’s 
sectional viscous drag is calculated based on either stripwise two-dimensional empirical drag 


coefficients or a stripwise two-dimensional integral boundary layer calculation. 


1.3. Euler/IBLT Throughflow Solver 


MTFLOW is a Euler/IBLT throughflow solver developed by Drela?. MTFLOW incorporates 
a streamline curvature method is a powerful method for solving the inviscid flowfield within a 
channel or annular passage. As the method is inviscid, a boundary layer solution must be used 
to calculate boundary layer quantities. The inviscid streamlines are then displaced off the body 


by the displacement thickness, é, which causes a redistribution of the inviscid streamlines. 


*Mark Drela, Associate Professor, Department of Aeronautics and Astronautics, MIT 
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Figure 1-1: Finite volume geometry for throughflow solver[2]. 


1.3.1 Euler Throughflow Solver 


MTFLOW is based on the conservative formulation of the steady state Euler equations[1]. It 
relies on a finite-volume throughflow method based on the streamline curvature method. Figure 
1-1 displays the finit-volume geometry representation used within MTFLOW. 

The first governing equation is the conservation of mass. Thus, across adjacent finite 


volumes boundaries 


P= —_ ox —_ 
PyUm181- Ai = poUme82- Ag (1.4) 


The conservation of momentum is written as the steady state Euler equation. 


— ~ => ~ a — ~ => - > 
P, A y+ (P:Um1 5] : Ay)Um181 +I BT = P» Ao + (pPoVm232 : A2)Um282 + It Bt (1.5) 


Conservation of energy is solved for in terms of enthalpy 


1 , 
A he he, (18) 


Additionally the average pressures along the faces pseudo-normal to and along the stream- 


tube must be equal. 


Pit P, =I +00 (1.7) 


MTFLOW further simplifies equations 1.4 through 1.7. The meridional flow speed, um, is 


obtained from a local streamtube conservation of mass. 


Q 


rr pA(2ar — BT¢) 


(1.8) 

The streamwise momentum equation solved by MTFLOW is of the same form as Equation 
1.4. 

dP + pUmdum + pugdvg + Pd(AS) — pd(AW) = 0 (1.9) 

The change in work, say due to the presence of a propeller, AW enters the formulation 


through the Euler turbine equation. 


AW = | Qd(rvg) (1.10) 


with 
; rad 
Q = blade rotation rate (—) 
sec 


This results in a computationally efficient methodology for solving the inviscid equations of 
motion, and is ideally suited to couple with the propeller lifting surface code PBD. 
1.3.2 Integral Boundary Layer Theory Equations 


Laminar and turbulent boundary layer flows are governed by the same integral momentum 
equation, obtained by integrating a combination of the continuity and simplified z-momentum 


equations term by term across the boundary layer [17]: 
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—+(24+ H)— = Lil 
d et ‘5, dx 2 aay 
where 
:  u u 
§ = momentum thickness = i —(1— —)dy (12) 
0 Us Ue 
* 2 U 
6 = displacement thickness = I (1—- ay (1.13) 
0 e 
6 
Hf = momentum shape factor = a (1.14) 
Ou 
Tw = wall shear stress = pv— |w (1.15) 
dy 
: ae : 2Tw 
C; = skin friction coefficient = a (1.16) 
€ 


In laminar flow, the variables 6, H, and C; can be reasonably well related with one-parameter 
velocity profile approximations. For marine vehicles moving at typical speeds, however, the 
transition from laminar to turbulent flow within the boundary occurs near the bow, and the 
boundary layer remains turbulent over the remainder of the vehicle. The turbulent-flow pro- 
file is complicated in shape, and many different correlations or additional relations have been 
proposed to effect closure of Equation 1.11 [17]. The closure relations used within MTFLOW 
are those derived by Swafford[15]. It is important to note again however, that MTFLOW does 
not solve for the boundary layer velocity profile, but rather computes the boundary layer pa- 
rameters and displaces the inviscid flow field off the body by 6 , and thus computes the global 
flow field inviscidly. 

One of the simplest, and thus robust, approaches to express the boundary layer velocity 


profile and recommended for general use by White [17], is the a power law: 


1/7 
cen ae 
Ue (i) ae 
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where ég9, the boundary layer thickness, is the value of y at which the boundary layer 
velocity is equal to 99 % of free stream velocity. By substituting Equation 1.17 into Equation 
1.13 and replacing the upper limit of infinity of the integral with d99, the relationship between 


6 and dgq is: 


6o9 = 86° (1.18) 


This relation is useful to determine the boundary layer thickess given a known displacement 


thickness. Then Equation 1.17 may be used to determine u as a function of y. 


1.3.3 MTFLOW Operation 


MTFLOW consists of three executable programs; MTSET, MTFLO, and MTSOL. MTSET 
creates the body geometry and automatically sets up the initial grid. MTFLO allows the 
input of external flow field information, such as rotor generated swirl or additional losses, into 
the flow domain. MTSOL executes the flow field solution. The only required input file for 
MTFLOW is the walls.xxx file which is used by MTSET to define the body geometry. The 
creation of the walls.xxx file and the execution of the MTFLOW programs are sufficiently 
covered in Reference [1] and will thus not be repeated here. Hereafter, MTFLOW refers to the 
entire throughflow solver, while MTSET, MTFLO, and MTSOL refer to the individual running 


of the respective executable programs. 


1.4 Objective 


This thesis provides a description of the fully automated coupling between the throughflow 
solver MTFLOW and the lifting surface propeller blade design and analysis code PBD. An 
overview of the process is first presented, followed by more detailed descriptions of the coupling 


processs. Finally, validation examples are presented. 
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Chapter 2 


Coupling PBD to MTFLOW 


2.1 Process Overview 


The coupling of the lifting surface propeller blade design and analysis code and the streamline 
curvature code follows the same logic as utilized by Kerwin et al [9] when coupling the propeller 
blade design and analysis code with a RANS code. The basis of Kerwin’s coupling was the use 
of distributed body forces within the RANS domain to represent the presence of the propeller. 
In the streamline curvature method, the propeller is modeled as a swirl introduction device. 
Thus added swirl, or angular momentum, is used to represent the presence of the propeller 
in the flowfield. This technique was first demonstrated by Renick [14] in a manual coupling 
method, and the present research is an extension of those efforts to produce a fully automated, 
robust method capable of solving the full range of propulsor types. While the overall concept 
and code execution is similar to that used by Renick, the coupling codes have been completely 
rewritten and additional changes have been incorporated into the flowfield solver to affect a 
correct and fully automatic coupling technique. 

The coupling technique is started by solving the nominal flowfield, z.e. the body with no 
propeller present. For viscous flow, the boundary layer is reconstructed and merged to the 
flowfield. This nominal flow field in then passed to the lifting surface code which solves the 
localized propeller problem. The updated circumferential mean induced velocities are converted 
to swirl, and are then inserted into the global flowfield which is then reconverged. The process 


is then repeated until both the flowfield, and lifting surface propeller results have converged. 
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PBD 
-Lifting Surface 


BL2Body 


-Boundary Layer 
Reconstruction 


~ PBDIMT 


-Input rV, and losses into 


a MTFLOW 


-Inviscid flowfield (to 5 *) 
-Boundary layer parameters 


Figure 2-1: Coupling sequence between propeller blade design code (PBD) and the throughflow 
solver (MTFLOW). 


The overall sequence is shown in Figure 2-1. 


2.2 Running the Coupling 


The overall coupling is controlled by the coupling administration file, mtcouple.inp. The 
coupling is executed by two computer codes: PBD2MT, which converts the PBD output 
to MTFLOW, and BL2BODY, which passes the updated flowfield to PBD. The codes are 
described in Table 2.1. 


The coupling administration file, mtcouple.inp, contains the required conversion informa- 


tion to execute the coupling. A sample administration file is displayed in Table 2.2 
The input lines have the following uses within the code: 
1. REYNOLD’S NUMBER: Reynold’s number for the throughflow domain. An inviscid 
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Description 


PBD2MT | 1) Converts the PBD output circumferential mean induced velocity to swirl. 
2) Writes out tflow.xxx, the input file to MTFLOW containing the 


induced rVp. 
BL2BODY | 1) Converts the MTFLOW output into velocity flowfield. 
2) For viscous cases, reconstructs the boundary layer profile. 


Table 2.1: Coupling code description 


Reynolds number 

inlet mach number 

Vship used in PBD to calculate Js 
x location of LE tip 

r location of LE tip 


MTFLOW case name 

pbd input file name 

BLinput toggle (O=no,1=yes) 

Nominal velocity toggle (0=no,1=yes) 
Number of blade rows 


OMAN AAP WHY eS 


Table 2.2: Example coupling administration file 


case is forced by specifing a Reynold’s number less than 10. 


2. INLET MACH NUMBER: For compressible flows, the Mach number is defined as the 
inlet velocity divided by the speed of sound. For water, M = zgghm. For general 
use with water a Mach number of 0.01 is recommended. This prevents round off errors 


internal to MTFLOW, and models the flow as virtually incompressible. 


3. VSHIP: This is the upstream velocity at infinity. For coupling purposes it is the fluid 
velocity at the domain inlet. The advance coefficient, J, within PBD is calculated using 


Vsnip- Traditionally cases are non-dimensionalized with Vnip = 1.0. 


4. X LOCATION OF LE TIP: The axial location of the propeller leading edge tip in 
MTFLOW coordinates. . 


5. R LOCATION OF LE TIP: The radial location of the propeller leading edge tip in 
MTFLOW coordinates. | 


6. MTFLOW CASENAME: The casename of the walls geometry file read by MTSET. 


The generated state and tflow.xxx file will be automatically labeled with the casename 
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10. 


extension. 


. PBD INPUT FILE NAME: Name of the PBD administration file. Used to read the 


advance coefhcient and number of blades. 


. BLINPUT TOGGLE: Specifies whether an inlet boundary layer profile is present. If set 


to 1, then the file BLin.txt is read, converted to entropy, and added to the flowfield. 


. NOMINAL VELOCITY TOGGLE: The toggle is set to 1 if the nominal profile is required 


for use in adjusting the specified input boundary layer to match a desired nominal profile 


at the plane of the propeller. 


NUMBER OF BLADE ROWS: Specifies the total number of rotor and/or stator rows 
present in the propulsor. Thus, a value of 1 describes a traditional propeller, and a value 
of 2 is correct for a pre-swirl/rotor combination. Up to a total of nine blade rows are 


permitted. 
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Chapter 3 


PBD to MTFLOW Conversions 


3.1 Program Overview 


Information is passed into MTFLOW by means of the tflow. xxx file. The tflow file may contain 
information for single blade rows, multiple blade rows, entropy loss due to the presence of an 
upstream boundary layer, or any combination of the above. This link of the coupling therefore 
requires the transfer of the necessary information into the tflow file. This is accomplished by 


the program PBD2MT. 


3.2 Program Operation 


3.2.1 Program Flow and Required Files 


After the successful execution of both MTFLOW to compute the nominal flowfield and the 
initial PBD run, PBD2MT is used to create the tflow.xxx file. As successive iterations are 
performed, grid information from the previous MTFLOW run, and the updated PBD output, 
are used to write an updated tflow file. By comparing the sizes of the various grid files, the 
program automatically determines the existence of an open, ducted, or internal flow propeller; 
and incorporates the necessary details into the tflow file. The required input files, source, 
and use, are explained in Table 3.1 The program operation is diagramed in Figure 3-1. The 
details of more complicated situations of multiple blade rows and inlet boundary layer profiles 


are described in later sections. 
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Uses propeller advance coefficient to compute rotation rate. 
mtcouple.inp Contains location of propeller in the flow domain and casename. 


PBDOUT.CMV Contains induced cimcumfential velocity. 
Gridxxxxx.dat | MTSOL | Contains grid information from previous iteration used to 
compute nondimensional streamline numbers. 


Table 3.1: Input files required by the program PBD2MT. 


a“ 
f, ie \ 
/ \ 
7 | MTSOL \| PBDI4 | 
/ ; User Created 
| GridLower.dat ] PBDOUT.CMV xxx.pbd 
\ GridUpper.dat / mtcouple.inp 
\ GridTotal.dat / 
ae | . / 
N Se el 4 Pa 
VA PBD2MT 
From previous iteration 
of MTFLOW 


tflow.xxx SWIRL.PBD 


Figure 3-1: Program order and file passing when running PBD to MTFLOW. The end product 
is the ascii file tflow.xxx, which is the required input for MTFLO. 


3.2.2 The tflow.xxx file 


Information is passed into MTFLOW by use of the tflow.xxx file which is read by MTFLO. 
Prescribed distributions of rVg, AS, etc., are specified as an MxN grid of points arranged in 
a logical rectangle in the s,t parameter space. In the current coupling routine, only rV¢ is 
specified in the file. The grid is input as a set of “profile blocks”. The grid is aligned with s 
the component perpendicular to the flow, and ¢ the eoatponent parallel with the flow. Each 
block has a constant s value, and thus defines quantities as a profile perpendicular to the flow 


at that s value [1]. An example tflow file is shown in Table 3.2. 


DH/a*2 


1.0 


Kstage(1) Nbilade(i) Omega(1) 
t(1,1) X(1,1) ¥(1,1) RVT(1,1) THK(1,1) DS(1,1) DH(1,1) 
t(2,1) X(2,1) Y(2,1) RVT(2,1) THK(2,1) DS(2,1) DH(2,1) 


t(N,1) X(N, 1) Y(N,1) RVT(N,1) THK(N,1) DS(N,1) DH(N,1) 


Kstage(M) Nblade(M) Omega(M) 


t(N,M) X(N,M) DAC, 


Y(N,M) RVT(N,M) THK (IN, M) Ds (i,m) 


Table 3.2: Example tflow.xxx input file to MTFLO 


The contents of the file are described below. Once again, items 11 and 12 are not determined 
in the current coupling technique, but are features which, in future work, may be utilized for 
further improvement. The reference length, L,e;, is defined as 1.0 in the MTFLOW coordinate 
system. ‘The points are specified with n values increasing perpendicular to the flow, and m 


values increasing parallel to the flow. 
1. LABEL: Arbitrary label text line. 
2. #...:. Comment label lines. 


3. *1.0...: Multiplier line which is applied to all subsequent values in the associated column. 


PBD2MT sets all multipliers at 1.0. - 


4. Kstage(m): Number of the rotor or stator stage which is described by the profile block. 
This MTFLOW function is not used in the coupling method, and is set to zero by 
PBD2MT. 


5. Nblade(m): The number of blades present in the current s line. 


6. Omega(m): The rotation rate in OLres/qint- 
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7. t(njm): The normalized streamfunction t, described below. 
8. X(n,m): The x coordinate defined as z/Lres. 
9. Y(n,m): The y coordinate defined as r/Lrez. 
10. RVT(n,m): The swirl defined as rVo/diniLref- 
11. THK(n,m): The blade thickness T9/Lref. 
12. DS(n,m): The added entropy AS from some adiabatic loss process. 


13. DH(n,m): The added enthalpy AH/a?,, from heat release. ( Not applicable to propeller 


inl 


analysis.) This is always set to zero by PBD2MT. 


3.2.3. Calculation of the Normalized Streamfunction 


The normalized streamfunction, t, is defined as: 


(3.1) 


where, in the MTFLOW application, w,and t, are the streamfunction values at the innermost 
and outermost streamsurfaces of the flowpath between any two bodies, or elements, in the 
flowfield. With multi-element configurations, y,and 2 are the streamfunction values at the 
bottom and top of one flow channel, and f is additionally offset by the channel number, so that 
0 <t <1 between the first and second element, 1 < t < 2 between the second and third element, 
etc. In all situations, each streamsurface has a unique t value associated with it [1]. For a 
given z,r location, the normalized streamfunction may be calculated by: 

p= taig PRT AE (3.2) 


Rmax 2 
f Row prr?dr 


which for uniform flow reduces to: 


R?- R. 
t= => OS (3.3) 
: ax Ae 


where 
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Z5, Re streamline n+1 
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Figure 3-2: Streamline tracing required to find nondimensional streamline number for any given 
point of interest in the flow field. 


yy 
Il 


radius of point of interest 
Rmin = radius of bottom of local flow field 


Rmax = radius of top of local flow field 


However, uniform flow generally only exits at the inlet of the MTFLOW domain. Thus, 
the following simple and robust approach is used to calculate the value of t given any z,r point 
of interest. The flowfield streamlines from the previous MTFLOW iteration are available in 
the Gridxxx.dat files. The applicable Gridxxx.dat file is selected, and the streamline above 
and below the point of interest is determined. The fractional difference between the point of 
interest and the identified streamlines, is traced back to the inlet, where uniform flow exists 
based on the inlet MTFLOW boundary conditions. Figure 3-2 shows the utilized streamline 


geometry. Thus the following equation is used to calculate R: 


6r1 


He 6r1 + dr Peas Bi (3.4) 
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Bz0 B=0 
Q+0 Q=0 


a _— ae 


Figure 3-3: Field parameter s,t grid for an open propeller on a straight shaft, showing doubled 
t profiles at blade leading and trailing edge. 


This calculated R value is then used in Equation 3.3 to determine the desired streamfunction, 


£, given any z,r position. 


3.2.4 Tflow.xxx Construction 


Figure 3-3 shows a general s,¢ grid for an open propeller on a straight shaft. To give more 
flexibility in splining, a slope break or discontinuity in either t or s can be specified by two 
identical t-value or s- value lines in each block of the tflow.xxx file. As shown in the example, 
the blade leading and trailing edge (s = 2 and s = 4) are doubled to force a spline break in 
B,Q, and rVg. Thus two identical profile blocks with identical X, Y are generated. One with 
B,Q, and rVg set to zero, and the other with B,], and rV@ set to the actual blade leading 
(or trailing) edge values. Similarly, a slope break or discontinuity in the s direction may be 
specified by two identical data lines within a given profile block. This is done at the blade tip, 


to enforce a zero rVg condition outside the propeller. 
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ee MTFLOW grid lines 4. _ PBD blade control points 


rV, input s lines @ Location of added rV, 


Figure 3-4: Sample grid demonstrating the relationship between the various parameters within 
PBD2MT. 


3.3. Multi-Blade Row Method 


In multiblade row applications, the induced rVg from one blade row is convected in the wake to 
the next blade row, at. which point the induced rVg from the next blade are added to those in 
the wake from the upstream blade. Kelvin’s theorem of conservation of circulation states that 
for any ideal fluid the circulation is constant along any streamline. Thus, lines of constant rVg 
are convected along each streamline. Since the MTFLOW grid conforms to the streamlines, 
this coupling method is ideally suited for multiblade row applications. Individual tflow.xxx 
files are created for each blade row in their own subdirectory as described in Section 3.2.4. Each 
tflow.xxx is then moved to the master case directory and renamed based on the order of the 
blade rows. For exmple, the first blade row tflow.xxx file is renamed tflow.B1, the second 
blade row tflow.xxx file is renamed tflow.B2, etc. The program BUILDTFLOW then takes 
all the tflow.Bx files, determines the cumulative rVg at each s line location, and creates the 
tfilow.casename file containing all the blade row information. It is then this tflow.xxx file 
which is read by the program MTFLO. The recommended directory structure and required file 
organization for a rotor/stator case is displayed in Figure 3-5. The program BUILDTFLOW 


can handle up to a total of nine blade rows. A script file for running a multi-blade row case is 
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ff 
Start Files 
mtcouple.inp 
walls.xxx 


1) Run MTSET 


: GridLower.dat | 
2) Run MTSOL : 


: GridTotal.dat i 


3) Run BL2BODY |} VELIOIN tec 


p> tHOW.B] rveersrsoreceesescerrenrsserans prvsvvssnennnnarasnscranesnnna tflow.B2 4, 
: Vv i 
10) Run BUILDTFLOW 

11) Run MTFLO 

12) Goto step (2) and repeat 


/ROTOR/ : H /STATOR/ 


Start Files i i Start Files 
mtcouple.inp i i mtcouple.inp 
rotor.pbd : i stator.pbd ; i 
rotor. BSN ¥ : stator.BSN v 
4) Run VELCON } _ | 7)Run VELCON 
5) Run PBD : : 8) Run PBD 


6) Run PBD2MT | 9)Run PBD2MT 


eee , Aenean : 
: Output 
tflow.stator { 


Figure 3-5: Multiple blade row coupling directory structure and program flow chart. The 
process is best controlled by a batch script file. 


contained in Appendix A. 

The final tflow.xxx file for a waterjet is shown in Figure 3-6. Due to the point density on 
the blades, the multiple points appear almost solid. Moving left to right through the flow field, 
the inlet s-line is shown first, followed by the closely spaced rotor s-lines containing the rotor 
generated rVp. The mid-stage s-lines represent the wake from the rotor. ‘The second set of 
closely spaced s-lines represent the stator, with their associated rV, being the superposition of 
the rotor wake and the negative swirl induced by the stator. The outlet section is the combined 
rotor/stator wake, and as the stator is designed to remove the induced swirl of the rotor, this 


final wake exhibits only a small amount of swirl. 
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Figure 3-6: Screen capture of MTFLO displaying the tflow.xxx file information for a ro- 
tor/stator waterjet. 


3.4 Upstream Boundary Layer 


An additional feature of the coupled PBD/MTFLOW method, is the ability to model only 
the stern portion of the body of interest. This has several advantages, some of which are 
helpful for rapid design trade-off studies, and some of which are required to successfully model 
vehicles with a significant portion of the propeller within a boundary layer. The idea is to 
only grid, and spend computational effort, on the stern of the vehicle. The MTFLOW domain 
inlet is located forward of the propeller by at least two propeller diameters. Experimental data 
indicates that the presense of the propeller does not affect the flowfield forward of this point. 
Thus this size flow domain is the minimum required to adequately capture the effective wake[5]. 
Once the approriate region is gridded by MTSET, the presence of a boundary layer at the inlet 
to the domain is modeled by adding an entropy loss in each streamtube. MTFLOW is then 
run in the inviscid mode, but with the boundary layer captured as an entropy loss. While 
the boundary layer is assumed to convect and evolve appropriately as the shape of the body 
changes, the additional viscous drag effects (which would normally thicken the boundary layer) 
occurring between the inlet plane and the propeller plane are not included in the inviscid mode. 
However, this effect is considered minimal in comparison to the overall accuracy and speed of 
this coupled method, especially given the fact that the action of the rotating blade rows will 


tend to accelerate the fluid and cancel much of the additional boundary layer thickness due to 
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friction. 

Once the nominal flow is calculated, the propulsor induced swirl may be added, and the 
total flow field/propeller problem solved iteratively. The speed of this “stern only” method 
permits a rapid analysis of multiple propulsor configurations for design space exploration. The 
other significant advantage of this feature, is that MTFLOW can now model the contribution to 
the flow field by the portion of the propeller within the boundary layer displacement thickness. 
If run as designed, a significant portion of the propulsor on a submerged body would be within 
the displacement thickness, which, by definition is outside the MTFLOW domain. The desired 
nominal profile may be provided by other flow solvers, such as RANS, experimental data, if 


available, or MTFLOW itself may be used to model the full body and run in the viscous mode. 


3.4.1 Velocity to Entropy Conversion 


To convert a velocity profile into an entropy distribution, we start with the entropy equation: 


(a) 
S=in|— 3.5 
Po is) 
For small changes (low Mach limit), 
spy. T (US a? 
Ree ont (Ut) (3.6) 
Po 2 Po 
Using the relation, 
pU? = ypM? (3.7) 


and, for the case of low Mach numbers (i.e. water) assume that the pressure is approximately 


equal to the absolute total pressure 


P*Po (3.8) 


so that 
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Iniet Boundary Layer Profile Specified as 
Entropy Loss on a Straight Shaft 
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Figure 3-7: Straight shaft with 4”” power law boundary layer specified by an entropy loss. 
& & 7 ?P y PY 


_ 1 (U2 -2?) 
AS =—57M ip (3.9) 


While the use of the specific heat ratio, y, and Mach number, M, with water flow are 
admittedly foreign, their use is important to correctly couple with MTFLOW. A specific heat 
ratio corresponding to that of air, y = 1.4, is required to ensure consistency with the built in 
MTFLOW methodology. The Mach number is as specified in the coupling administration file, 
mtcouple.inp. 

To validate the use of entropy as a boundary layer modeling tool, a straight shaft test case 
was conducted with a ay power law boundary layer present. The results are demonstrated in 
Figure 3-7. The agreement between the desired boundary layer, and that convected through 
MTFLOW, displays the effectiveness of simulating a boundary layer as an entropy loss. 
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3.4.2 Boundary Layer Entropy Modeling on Submerged Body 


When the same procedure is conducted on a body of revolution, the boundary layer profile at 
the domain inlet matches the specified profile, but the boundary layer downstream develops dif- 
ferently than the experimentally measured profile. This is presumably a result of the increasing 
invalidity of.the pressure assumption of Equation 3.8 as the curvature of the body increases. 
However, to ensure that only velocity data is required for the case setup, it is not desirable 
to introduce a pressure correction. ‘To overcome this deficiency, a procedure is available to 
tune the specified inlet velocity profile, and the associated added entropy loss, to result in the 
correct nominal flow field in the propeller plane. Thus, while the resultant inlet flow domain 
may differ from the measured profile there, as the flow develops the correct nominal profile is 
present in the region of the propeller. Once the necessary inlet profile is determined, it remains 
constant for use with the propeller present. 

The procedure and program used to scale the inlet boundary layer to achieve the desired 
nominal profile at the propeller is displayed in Figure 3-8. Haung Body 1 (DTMB Model 
5225-1) was used to validate this approach [8]. The file format of the boundary layer profiles, 
BLin.txt at the domain inlet and Exp977nom.dat at the propeller plane, are shown in Table 
3.3. MTFLOW will correctly determine the radial velocity component based on the body 
shape, so only the axial velocity component must be specified in the boundary layer profile 
files. The boundary layer AS loss, is added in the same manner as the blade induced rVg and 
appears as an MTFLO input variable DS as shown in Table 3.2. The BLINPUT TOGGLE in 
the coupling administration file mtcouple.inp must be set to 1 to include the AS loss. During 
the setting of the inlet profile to match the desired nominal profile at the propeller plane, a 
PBDOUT.CMV file is required to specify the blade location; however, no rV¢ is desired. This too 
is controlled by the mtcouple.inp file by setting the NOMINAL VELOCITY TOGGLE to 1 
during this scaling. To prevent overshooting the desired profile, a damping coefficient of 0.25 


is included in each scaling step. The scaling is accomplished for each streamline as follows: 


1 /V; Vets 
Vi iia f _i (Anita (3.10) 
4 VurrFriow () 
where 
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V(1) 
v(2) 
v(3) 


V(NumPts) ) 


Table 3.3: Required format of boundary layer file. The number of points in the boundary layer 
are listed in the first line, followed by each R, V point, where V is the axial velocity component. 


Varin(Y)new = updated streamline axial velocity in boundary layer input file BLin.txt 
Verin(~)od = previous streamline axial velocity in boundary layer input file BLin.txt 
Vurriow(?) = streamline axial velocity at propeller plane in MTFLOW 
Voesirea(¥) = Secale desired axial velocity at propeller plane in file Exp977 .nom 


Validation of this method was conducted with Huang Body 1 (DTMB Model 5225-1). Fig- 
ure 3-10 shows the required inlet profile which results in the current boundary layer in the 
sroweller plane shown in Figure 3-9. The complete stern section of the body and boundary 
layer growth is shown in Figure 3-11. Of note, a far field setting of a constant pressure jet 
boundary in MTSOL was required for stable convergence. This is imposed by setting the far 


field type to “3” in the solution parameter settings. 


3.4.3 Gridding Requirements for Stern Section Modeling 


The automatic gridding procedure contained in MTSET was designed for wings, ducts, or 
internal flow. Consequently, when directly applied to the stern section of a submerged body, 
the grid smoothing routines may have difficultly in the inlet region adjacent to the body. The 
bottom left corner may drift above the next streamline as shown in the top of Figure 3-12. By 
introducing a spline break as a doubly specified body point as demonstrated in Reference [1], 
and artificially lowering the inlet body point slightly, correct gridding is accomplished. This 
ensures the bottom streamtube has positive area, and thus an initial positive velocity. The 


initial grid is viewed in the TECPLOT formatted file ORGGRID.tec. The existence of a negative 
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mtcouple.inp 
PBDOUT.CMV * 


Entered Profile 


BLinax PBD2MT 
MTFLO/MTSOL/BL2BODY 


VELJOIN.tec 


MT977nom.dat ** 


! Desired Profile 
ScaleBL  |«—————Exp977nom.dat 


| eee er Oe ee ee ee BLin.txt * Uses blade location data only 


** Created using TecPlot polyline extraction 


Figure 3-8: Program operation and file passing to scale the specified entropy loss to match the 
desired nominal profile at the propeller inlet plane. 
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Nominal Velocity Comparision 
Huang Body 1 at X/L= 0.977 
Inlet Entropy Set for Desired X/L = 0.977 Profile 
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Figure 3-9: MTFLOW and experimental nominal velocity comparison of Huang Body 1 at X/L 
= 0.977 with inlet entropy loss specified at X/L = 0.914. 
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Nominal Velocity Comparision 
Huang Body 1 at X/L= 0.914 
Inlet Entropy Set for Desired X/L = 0.977 Profile 
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Figure 3-10: Resulting nominal velocity profile at X/L = 0.914 in MTFLOW with the required 
entropy loss to correctly match the nominal profile at the propeller location of X/L = 0.977 


compared to experimental profile at X/L = 0.914. 
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Figure 3-11: Nominal flow over stern section of Huang Body 1 with inlet boundary layer modeled 
as an entropy loss. 
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Inlet Body Points Entered at Actual Location 


20.06 ~ 20.07 a 20.09 20.1 


0.06 20.08 : 70.08 204 


Figure 3-12: Grid comparsion showing drift of bottom left corner (top) when inlet body points 
entered at actual location versus adjusted body points (bottom). 


area can at a minimum result in localized flow difficultly, and in the extreme, can lead to a 
negative temperature and the subsequent crashing of the program MTSOL. Fortunately, once 
a successful grid is created, it is sufficient for all subsequent coupled runs. A corrected grid is 
shown in the bottom half of Figure 3-12. The localized point adjustment has no impact to the 


overall flow downstream in the propeller region. 
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Chapter 4 


MTFLOW to PBD Conversions 


4.1 Program Overview 


The link from MTFLOW to PBD is accomplished by the program BL2BODY. The output 
from BL2BODY is the Tecplot® formatted flow field velocity file, VELJOIN.tec, which is 
read by VELCON [4], the traditional velocity conversion program used to create the updated 
flow field for use by PBD. For an inviscid case, BL2BODY simply converts the MTFLOW 
output to the correct format required by VELCON. In a viscid case, the MTFLOW output 
only contains inviscid velocities for a body offset by the displacement thickness, and thus the 
flowfield must be expanded to the body boundary, and an appropriate boundary layer velocity 


profile reconstructed. 


4.2 Program Operation 


The overall program execution order is displayed in Figure 4-1. The coupling administration 
file, mtcouple.inp, is read to determine whether a viscid or inviscid case exits. If the specified 
Reynolds number is less than ten, an inviscid case is assumed. For the inviscid case, the 
VELJOIN.tec file is written with the velocities referenced to free stream velocity, based on the 
velocity and Mach number as referenced in the coupling administration file. For a viscid case, 
a boundary layer reconstruction occurs. The beginning and ending points of the bodies, and 


the number of total bodies in the flow field, are determined from the walls.xxx file. The body 
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Figure 4-1: Program order and file passing when running MTFLOW to PBD. The end product 
is the TECPLOT format file restart.vel (or other name as specified during VELCON execution), 
which is the required input for PBD. 


geometry is read from the ORRIGGRID.tec file. As this file is created by MTSET, before a flow 
solution is created, this grid conforms to the actual body location. The original work by Renick 
[14] utilized a Swafford boundary layer profile, however, the present coupling reverted to a pure 
a power law reconstruction. While this is less accurate in high pressure gradient regions, 
it is more robust, and was thus more consistent with the goal of a rapid, stable, automatic 
coupling method. Certainly, the addition of a more advanced, yet robust boundary layer 
reconstruction technique would improve the coupling accuracy. Of note, the recommended, and 
often necessary method of characterizing a boundary layer due to the presence of a significant 
displacement thickness in the propeller region for a viscid case, is the entropy loss method 
described in Section 3.4.2. When a boundary layer is created using an entropy loss, MTFLOW 
is run in the inviscid mode, and the drawbacks of the 1 power law are irrelevant. 


The velocity output from MTSOL is written in the file OUTVEL.tec, and for a viscous case, 


the boundary layer information is written in the file OUTBL.tec. A sample OUTVEL.tec velocity 
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field is shown in Figure 4-2. The maximum displacement thickness, maz 6*, is determined from 
the OUTBL.tec file, and following the derivation in Section 1.3.2, the maximum boundary layer 
thickness, maz 699, is calculated by Equation 1.18. The streamline immediately outside the 
maz 699 is considered the lowest fully inviscid streamline. The MTFLOW grid is cut at this 
location, and boundary layer reconstruction is conducted between the body and the lowest fully 
inviscid streamline. First, new high density grid lines are faired in between the lowest fully 
inviscid streamline and the body. Then, velocities are assigned to the new grid nodes. If an 
individual node is above the local 699, then the OUTVEL.tec grid is interrogated to determine 
the local velocities. If an individual node is between the body and the local 699 then the local 
axial and radial velocity components are calculated in accordance with the ym power law of 
Equation 1.17. Between the local 6* and local é99 the tangential velocity is assigned the value 
in the MTFLOW domain as contained in OUTVEL.tec. Between the body and the local 6 
no tangential velocity data exists within MTFLOW. To resonable approximate the tangential 
velocity in this region, the tangential velocity is linearly interpolated between the actual value 
at the edge of the MTFLOW domain, 6*, and an assumed value at the wall of 50 % of that at 
the edge. Figure 4-3 shows an example VELJOIN.tec file written by BL2BODY with the new 
viscid, dense grid joined to the MTSOL calculated fully inviscid grid region. 

To compare the boundary layer reconstruction with the experimental data, velocity profiles 
were extracted from the converged viscous Huang Body 1 nominal profile case. Figure 4-4 
shows the results at four X/L locations. As expected, the TP power law fails to accurately 
represent the experimental profile in the tapered stern section. However, it still serves as a 


useful tool in the absence of more detailed information, and for use along straight shafts, ducts, 


and internal flow cases. 
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Figure 4-2: Sample viscous flowfield output from MTFLOW. 


= faired 


Figure 4-3: Sample BL2BODY grid with boundary layer reconstruction grid added to the 
inviscid MTFLOW grid. 
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Figure 4-4: Comparison of reconstructed boundary layer profile and experimental measurements 
for Huang Body 1. As expected, the accuracy of the a power law diminishes in the highly 
tapered stern region. 
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Chapter 5 


Validation 


Code validation was conducted on the following cases: 


1. Open Propeller 
2. Ducted Propeller 
3. Waterjet 


4. Submerged Body 


The open water comparison used propeller 4119 in comparison with the 1998 International 
Towing Tank Committee (ITTC) tests[13]. The ducted propeller KA-455 from the Netherland 
Ship Model Basin (NSMB), Kaplan series ducted propulsor [16], was utilized for the ducted 
case. Internal flow was validated by the WaterJet-21, as tested in the Marine Hydrodynamics 
Laboratory (MHL) at the Massachusetts Institute of Technology (MIT){10]._ Huang Body 1 
(DTMB Model 5225-1) with propeller 4577 [5] was used for the submerged axisymmetric body 


comparison. 


5.1 Open Propeller 


Propeller 4119 on a straight shaft was used as the initial validation of the coupling method. 
Additionally, the final results were compared with the experimental results from the 1998 ITTC 


tests as the open propeller test case [13]. The inviscid axial effective velocity was calculated 
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and is show in Figure 5-1. Given a nominal axial inflow of 1.0, the axial effective inflow 
velocity should be equal to 1.0 on the entire blade. The near unity results demonstrate the 
effectiveness of the coupling method. An additional useful coupling check is to compare the 
tangential trailing edge velocities between PBD and MTFLOW. This is synonymous with 
comparing the circulation distributions. Thus, the trailing edge circumferential mean induced 
tangential velocities output from PBD are compared to the tangential velocity at the same 
spacial location in the MTFLOW galeilated flow domain. This is shown in Figure 5-2. The 
exact agreement between both domains is an important verification that the coupling is being 
carried out correctly. Propeller 4119 results are compared with the 1998 ITTC experiment 
results in Figure 5-3. MTFLOW was run in the viscous mode, allowing boundary layer growth 
on the shaft. PBD was run with both the sectional drag coefficients and thickness effects 
included. A relatively dense lifting-surface grid of 35 x 35 vortices was required for grid 
convergence!. As expected, the accuracy of the solution diverges at off design J values due to 


the inability of the lifting surface method to adequately handle extreme blade angles-of-attack. 


5.2 Ducted Propeller 


A ducted propulsor validation was conducted using the Kaplan KA-455 with nozzle 20 as tested 
in the Netherlands Ship Model Basin{16]. For comparative purposes, the thrust value contains 
only the rotor generated thrust, and does not contain the nozzle effect. Figure 5-4 is the 
KA-455 as modeled in the lifting surface tool. MTFLOW was run for both an inviscid and 
viscid case. Once again, the PBD analysis included both the sectional drag coefficients and 
thickness effects. The results are shown in Table 5.1. The results demonstrate the accuracy 
of the method and the minimal effect that the viscous boundary layer on the shaft and duct 


have on the overall solution. 


1Grid convergence means that a further increase in the grid density does not affect the results. 


AT 


P4119 Inviscid Axial Effective Velocity 


Figure 5-1: Countor plot of the effective axial velocity in the presense of a uniform nominal 
inflow velocity of 1.0 for P4119 on an inviscid straight shaft. 


4119 Trailing Edge Vtheta Comparison between 
PBD and Mtflow 


——t}-— PBD Swirl 
——*— MTFLOW Swirl 


Figure 5-2: Comparison of the circumferential mean tangential velocity at the blade trailing 
edge of P4119 from MTFLOW and PBD. 
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4119 Coupled MTFLOW vs 1998 ITTC Experiment 
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Figure 5-3: Comparison of P4119 MTFLOW/PBD results with 1998 ITTC Experiment. PBD 
was run with both sectional drag coefficients and thickness effects included and with a vorticy 


grid lattice of 35 x 35. 


Figure 5-4: Kaplan KA-455 ducted propulsor as modeled in PBD. 
MTFLOW/PBD | MTFLOW/PBD 


<a 
Re | invicid [500000 =| 


Table 5.1: Comparison of KA-455 inviscid and viscid MTFLOW/PBD rotor results with ex- 


perimental data at design J=0.36. 
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Figure 5-5: Contour plot of circumfential mean tangential velocity of WaterJet WJ21 with 
streamlines superimposed. Flow is from left to right. The rotor leading edge is located at 
approximatly X = 0.2 and the stator leading edge is located at approximatly X = 1.0. 


5.3 WaterJet 


Internal flow verification utilized the Waterjet WJ21 tested in the Marine Hydrodynamics 
Laboratory at the Massachusetts Institute of Technology[10]. The torque and thrust of the 
rotor were compared to an equivalent RANS calculation and with the experimentally measured 
torque. As with the ducted case, the thrust of the rotor does not correspond to the net thrust 
produced by the propulsion system, but does serve as a useful quantitative comparison. A plot 
of the circumferential mean tangential velocities with streamlines superimposed is displayed in 
Figure 5-5. The results, including a computational time comparison, are shown in Table 5.2. 
This particular case demonstrates the relative usefulness of the MTFLOW/PBD method. The 
results are quite close, in spite of the fact that this is a particularly complex flow case. The 
reduction in computing time from 12 hours to 30 minutes is very substantial. As the flow in 
the test section of the waterjet is dominated by potential flow, this case was run in the inviscid 
mode. For design applications, the boundary layer growth of the upstream flowpath would be 
perfectly suited to be modeled via an entropy loss. This would permit the blade design to be 
conducted within the actual flow field expected, while not requiring the complete gridding of 


the piping system as is currently required by RANS methods. 
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[| MTFLOW/PBD [RANS/PBD | Experimental | 
[Kr [192 CT ae TCid 
0.409 0.396 0.424 
12 hours 


Time 30 min 


Table 5.2: Comparison of waterjet WJ21 MTFLOW/PBD and RANS/PBD results with ex- 
perimental data. 
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Figure 5-6: Representation of Huang Body 1 (DTMB Model 5225-1). 
5.4 Submerged Body 


This test case is one of a series of axisymmetric bodies tested by the David Taylor Model Basin 
by Huang et al [7], [8], [6], [5]. The same forebody (DTMB Model 5225) was used for all 
experiments while a series of afterbodies with an increasing degree of taper. These bodies have 
been used extensively for validating solutions to the effective wake problem using analytic and 
numerical methods. The afterbody considered here, Afterbody 1, is a non-separating stern 
with a low tailcone angle. A profile view of the experimental body (DTMB Model 5225-1) is 
shown is Figure 5-6. It was tested in the presence of an open rotor in wind tunnel and towing 
tank facilities. An existing seven-bladed propeller (P4577), with a diameter of 54.5 % of the 
hull diameter is mounted at z/L = 0.983. As the actual body boundary layer displacement 
thickness includes a significant portion of the propeller, the entropy loss method of Section 3.4.2 
is used in this application. 

To correctly set the nominal profile, experimental boundary layer profiles were selected at 
two locations. The first, located at X/L = 0.914, served as the initial inlet plane profile, and 
the second, located at the propeller inlet of X/L = 0.977, served at the goal nominal profile. 
The inlet plane profile was then adjusted in accordance with the procedure of Section 3.4.2 until 


the MTFLOW nominal propeller plane inlet profile matched the experimental profile. Once set, 
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Nominal and Total Velocity Comparison of 
MTFLOW and Experiment 
Huang Body 1 with P4577, J= 1.25, at X/L = 0.977 
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Figure 5-7: Huang Body 1 nominal and effective velocity comparison of MTFLOW and exper- 
iment at X/L = 0.977 with inlet entropy loss specified at X/L = 0.914. The propeller leading 
edge tip is located at R = 0.545. 


the inlet profile remained fixed during the propeller analysis. The nominal and total propeller 
plane profiles are shown in Figure 5-7. The calculated and experimental thrust and torque 
coefficients are contained in Table 5.3. Of note, the listed experimental torque coefficient is 
not purely a measured value. The propulsor used in the self-propulsion experiment was a 
stock propulsor and was not operated at its design angle of attack. The experimental Kg was 
calculated by a propulsor performance prediction computer program based on hydrodynamic 
pitch angles tan @; which were iteratively scaled based on the ratio of measured to predicted 


Kr(6). 
A propulsor performance prediction computer program was used to compute the values of non-dimensional 
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Kr Kg | % Error | % Error 
Kr Ko 


Experimental 0.2276 | 0.0453 
MTFLOW /PBD-14.3 | 0.2392 | 0.0436 5.1% -3.7% 


Table 5.3: Comparison of propeller 4577 performance on Huang Body 1 with J = 1.25. 


The level of accuracy for Huang Body 1 is consistent with the use of PBD/MTFLOW as a 
rapid, early-stage, design space exploration tool. Capturing the effective wake problem in this 
manner allows the rapid analysis or design of multiple propellers in minutes versus hours, as is 


expected with RANS methods. 


circulation G and hydrodynamic pitch angle tan §, for the estimated values of ue. The final values of tan G; were 
scaled up or down by the ratio of the measured value of Kr to the computed value of Kr. The values of G and 
modified tan 8, were then used to compute the propulsor-induced velocities by a field-point velocity computer 
program with a lifting-surface option. The new values of u, were then used to compute a second estimate of 
the effective velocity ue. Three iterations were sufficient for convergence(4]. 
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Chapter 6 


Conclusions 


6.1 Summary of Results 


The coupled PBD/MTFLOW propeller blade design and analysis method is a viable alterna- 
tive to axisymmetric RANS codes for solving circumferential mean throughflow problems at a 
significantly reduced computational, as well as user-preparation, time. The limited number of 
required starting files and the automatic gridding routine of MTSET permits rapid problem 
set-up. The coupling works for open, ducted, and internal flow cases; for both single and 


multiblade rows. 


6.2 Recommendations for Future Work 


1. Addition of blade blockage effects in MTFLOW. A feature of the thlow.xxx file allows 
blade thickness to be passed into MTFLOW to model blockage, and thus this is an ideal 
feature to add to the PBD2MT program. The blockage effect becomes important for 


internal multi-blade propulsors. 


2. Addition of viscous blade drag in MTFLOW. In the current coupling method, the vis- 
cous blade effects are accounted for by PBD in the torque and thurst calculations, however 
the effect of this loss on the fluid velocity is not passed back into MTFLOW. The effect 
of the loss is small on single-blade row open propellers, but becomes significant for multi- 


blade row internal flow propulsors. Accounting for this loss would improve the accuarcy 
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of the solution in general, and would permit a realistic design capability for downstream 


blade rows in particular. 


. Fully automated inlet boundary layer profile modeling. The current inlet boundary layer 


profile method requires some user iteration to establish the entropy loss necessary to 
achieve the desired nominal profile. Altering the current scaling technique to a fully 
automatic version, transparent to the user, would further reduce the required user set-up 
time and level of difficulty. This is particularly beneficial in that the modeling of the 
stern section of a vehicle is presumably one of the most useful aspects of this coupling 


method, and is thus a feature which will presumably be frequently utilized. 


. Improved gridding for stand alone stern section. The gridding routine requires some 


manual iteration as discussed in Section 3.4.3. Eliminating this necessity would reduce 


the required set up time for stern section modeling. 


. Coupled Propeller Lifting Line (PLL)/MTFLOW. Coupling PLL and MTFLOW would 


result in a powerful design tool. The optimum circulation distribution could be designed 
by PLL in the effective wake as calculated by MTFLOW. Then, the blade shape design 
could be accomplished with PBD/MTFLOW. 
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Appendix A 


Multi-Blade Row Script File 


Examples 


This appendix contains example script files for a multi-blade row case. A script file for a single 
blade row application is a simplified version of the BatchFile below and is thus not included 


here. 


A.1 BatchFile 


TESTDIR=$HOME/RUNS/WJ 
ROTOR=$TESTDIR/ROTOR 
STATOR=$TESTDIR/STATOR 

bl 2body=$HOME/CODES/MTCouple/bl2body 
velcon9=$HOME/CODES/VELCON/VELCON9/velcon9 
pbd143=/nfs/propnut /CODES/CODES/PBD14/PBD14 .3/SRC/pbd14 
pbd2mt=$HOME/CODES/MTCouple/pbd2mt 
mtflo=$HOME/CODES/Mtflow/bin/mtflo 
mtsol=$HOME/CODES/Mtflow/bin/mtsol 
buildtflow=$HOME/CODES/MTCouple/buildtflow 
cd $TESTDIR 


rm Batch.log 
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date > $TESTDIR/Batch.1log 

pwd >> $TESTDIR/Batch.log 

MAX=10 

COUNT=1 

echo ’TITLE = ??WJ21 PBD/MTFLOW STATOR COUPLING’’’ >> $TESTDIR/CTQOUT.TOT 
echo °VARIABLES = ?’N’’, °?Ct’?’?, ??10Cq’’? >> $TESTDIR/CTQOUT.TOT 
echo ’?ZONE T=’’?Ct / 10Cq’’, I= ? $MAX >> $TESTDIR/CTQOUT.TOT 

echo ’TITLE = ??WJ21 PBD/MIFLOW ROTOR COUPLING’’’ >> $TESTDIR/KTQOUT.TOT 
echo ’VARIABLES = ??N??, ??Kt’’?, ??10Kq’?’ >> $TESTDIR/KTQOUT.TOT 
echo ZONE T=?’?Kt / 10Kq’?, I= ? $MAX >> $TESTDIR/KTQOUT.TOT 

cd $TESTDIR 
HHHHHHHHHAHHHHHRSEHRAEERERERERHEHRAAAHBR HHH ERAHR RRR ARRAS HR AHH 
while test $COUNT -le $MAX 

do 

echo ’ ## PBD - MIFLOW ITERATION ” $COUNT >> Batch.log 

echo ’ ####### PBD - MTFLOW ITERATION ° $COUNT’ #######? 
HEHHRHHHAPH HEHE HERR HHH ERHARRHRAR RHR RAHA R HARRAH HERR H RR HERR RE HE 

# MTFLOW TO PBD # 
HHFHHHHHHEHHAHHAEHRHRA EEA HEHHH RH RRHH APH RHA AHHHR AH RHEE RHA BHR ERAS 

# Create VELJOIN.tec from Mtflow results 

$bl2body >> Batch.log 

echo ’ ## Done with bl2body ##’ 

# Run velcon and pass req’d files to ROTOR and STATOR directories 
$velcon9 < velcon.rot >> Batch.log 

mv rotor.vel $ROTOR/. 

cp GridLower.dat $ROTOR/. 

cp GridUpper.dat $ROTOR/. 

cp GridTotal.dat $ROTOR/. 

cp VELJOIN.tec $ROTOR/. 


$velcon9 < velcon.stat >> Batch.log 


oT 


mv stator.vel $STATOR/. 

cp GridLower.dat $STATOR/. 

cp GridUpper.dat $STATOR/. 

cp GridTotal.dat $STATOR/. 

cp VELJOIN.tec $STATOR/. 
HHHHHHHHHHHHHHAHHHHHHH REHASH EHH R ARETE HERR ARR ERR RH AER R RES 
# PBD / PBD to MTFLOW # 
HHHHHHHHHHHHHHHHAHHHHHREHHHEEREREHHER EERE HASH RAB ES RRR ER RRR H 
# Run PBDi4.3 on Rotor and create tflow.rotor 

cd $ROTOR 

$pbd143 < pbd.in >> ../Batch.log 

echo ’ ## PBD complete on rotor ##’ 

$pbd2mt >> ../Batch.log 

cp tflow.rotor $TESTDIR/tflow.B1 

#cat PBDOUT.CMV >> PBDTOT.CMV 

#cat PBDOUT.CMF >> PBDTOT.CMF 

#cat PBDOUT.SGR >> PBDTOT.SGR 

cat PBDOUT.KTQ >> ktrot.tot 

tail -1 PBDOUT.KTQ > $TESTDIR/ktq.in 

read word1 word2 word3 word4 < $TESTDIR/ktq.in 
echo $COUNT $word3 $word4 >> $TESTDIR/KTQOUT.TOT 
# Run PBD14.3 on Stator and create tflow.stator 
cd $TESTDIR/STATOR 

$pbd143 < pbd.in >> ../Batch.log 

echo ’ ## PBD complete on stator ##’ 

$pbd2mt >> ../Batch.log 

cp tflow.stator $TESTDIR/tflow.B2 

#cat PBDOUT.CMV >> PBDTOT.CMV 

#cat PBDOUT.CMF >> PBDTOT.CMF 

#cat PBDOUT.SGR >> PBDTOT.SGR 
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cat PBDOUT.KTQ >> ctstat.tot 

tail -1 PBDOUT.KTQ > $TESTDIR/ktq.in 

read word! word2 word3 word4 < $TESTDIR/ktq.in 

echo $COUNT $word3 $word4 >> $TESTDIR/CTQOUT.TOT 
HEPAHHHEHHAHHEHRAAERH AHHH H REAPER HEHE REE EHH HHP AREA HR ES 
# MTFLOW # 
HHHHRHHHHEHHABHHEHHBEERERREPEH EERE ERR ER EE RRR R HH BH HH 
# Combine tflow files and run Mtflow 

cd $TESTDIR 

$buildtflow >> Batch.log 

echo ’ ## wrote tflow.wj21 ##° 

$mtflo wj2lexps < runMTIFLO >> Batch.log 

$mtsol wj2lexps < runMTSOL >> Batch.log 
HEHHHHHHHHHHHHHERHHHHHHRHH HEE HERERREBHPEE RBH ERE R EERE ARERR EH 
COUNT=‘expr $COUNT + 1° 

# End of while loop 


done 
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A.2 RunMTFLO 


P 


r 
{return] 
{return] 
Ww 


q 


[return] 


A.3 RunMTSOL 


WwW 


q 


[return] 


A.4_ velcon.rot 


VELJOIN.tec 
rotor.vel 
0.12396 0.52176 
1 

10 
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